Many-body effects in nuclear structure 
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We calculate, for the first time, the state-dependent pairing gap of a finite nucleus ( 120 Sn) diag- 
onalizing the bare nucleon-nucleon potential (Argonne 1114) in a Hartree-Fock basis (with effective 
A;— mass nik ~ 0.7 m), within the framework of the BCS approximation including scattering states 
up to 800 MeV above the Fermi energy to achieve convergence. The resulting gap accounts for 
about half of the experimental gap. We find that a consistent description of the low-energy nuclear 
spectrum requires, aside from the bare nucleon-nucleon interaction, not only the dressing of single- 
particle motion through the coupling to the nuclear surface, to give the right density of levels close 
to the Fermi energy (and thus an effective mass m* m m), but also the renormalization of collective 
vibrational modes through vertex and self-energy processes, processes which are also found to play 
an essential role in the pairing channel, leading to a long range, state dependent component of 
the pairing interaction. The combined effect of the bare nucleon-nucleon potential and of the in- 
duced pairing interaction arising from the exchange of low-lying surface vibrations between nucleons 
moving in time reversal states close to the Fermi energy accounts for the experimental gap. 



In the study of finite many-body systems such as the 
atomic nucleus with its rich variety of quantal size ef- 
fects, structural properties, and fluctuations, the central 
problem has been to identify the appropriate degrees 
of freedom for describing the phenomena encountered. 
The complementary concepts referring to the indepen- 
dent motion of the individual nucleons and the collective 
behaviour of the nucleus as a whole provide the elemen- 
tary modes of excitation needed to describe the system 
lj. The unifying picture emerging from the interweav- 
ing of these degrees of freedom and described in terms 
of nuclear field theory (NFT) [2-6] based on the particle- 
vibration coupling (for other related approaches, see e.g. 
[7-10]) and tailored upon QED [11-12], has been applied 
to a number of schematic models and realistic situations 
[13-21] and its validity demonstrated. It thus provides a 
natural framework to assess the role different degrees of 
freedom play in the nuclear structure. 

An important subject presently under intensive study 
concerns the characterization of an eventual long range 
component of the pairing interaction in nuclei [22-24]. In 
what follows we use NFT to assess the importance the 
exchange of vibrations between pairs of nucleons moving 
in time reversal states have in building up pairing cor- 
relations in nuclei, taking also into account self-energy 
and vertex corrections (i.e. avoiding using approxima- 
tions which, in condensed matter literature are connected 
with the so-called Migdal theorem, cf. ref. [25] and refs. 
therein) . 

To this scope, we study the quasiparticle and vibra- 
tional spectrum of odd- and even- isotopes of single- 
closcd-shell nuclei, where all the richness of the single- 
particle and collective degrees of freedom are fully ex- 
pressed, avoiding the extra complications of static de- 
formations and associated rotations. The spectra of 



the ^ Sn isotopes, in particular those with mass num- 
ber ^4=119,120 and 121, with their abundance of de- 
tailed experimental information, provide an excellent lab- 
oratory where to test the importance of the residual 
pairing interaction and its relation to self-energy pro- 
cesses. To be remembered that BCS theory connects the 
mass enhancement factor A associated with the uj— mass 
(m u = m(l + A), cf, e.g. [26]) with the pairing gap 
A 2uj c exp(— 1/A), oj c being the energy associated with 
a typical vibration (boson) of the system [25]. 

In general one fixes the parameters of the effective in- 
teraction of nucleons in the nucleus, by requiring mean 
field theory, as a rule Hartree-Fock or Hartree-Fock- 
Bogoliubov theory if the system is superfluid, to repro- 
duce the experimental findings: binding energies, mean 
square radii, etc. This is equivalent to requiring that 
the solution of the Schrodinger equation describing the 
bound states of the electron-proton system, interacting 
through the Coulomb force, reproduces the energy levels 
of the hydrogen atom. We know that this is not pos- 
sible, unless the renormalization effects arising from the 
electron-photon coupling are properly taken into account 
as prescribed by QED [11,12]. Similarly, the parameters 
of the effective nuclear interaction should reproduce the 
experimental findings only when the particle-vibration 
coupling is allowed to renormalize, screen and dress the 
different modes of elementary excitation and the inter- 
action among them in a similar way as, for example, 
the Lamb shift of the hydrogen atom is accounted for 
only when the renormalization effects arising from the 
electron-photon coupling are considered. 

The formalism we shall use is based on the Dyson equa- 
tion [22]. It can describe on equal footing the dressed 
one-particle state 5, of an odd nucleon renormalized by 
the (collective) response of all the other nucleons (Figs. 
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l(a)-l(d)), the renormalization of the energy fiuj v (Figs. 
2(a)-2(b)) and of the transition probability B{E\) (Figs. 
2(c)-2(f)) of the collective vibrations of the even system 
where the number of nucleons remains constant (corre- 
lated particle-hole excitations), and the induced interac- 
tion due to the exchange of collective vibrations between 
pairs of nucleons [23] , moving in time reversal states close 
to the Fermi energy (Figs. l(e)-l(g)). We include both 
self-energy and vertex correction processes, thus satisying 
Ward identities (cf., e.g., [25]). Within this framework, 
the self-consistency existing between the dynamical de- 
formations of the density and of the potential sustained 
by "screened" particle-vibrations coupling vertices leads 
to renormalization effects which make finite (stabilize) 
the collectivity and the self-interaction of the elementary 
modes of nuclear excitation, in particular of the low-lying 
surface vibrational modes, providing an accurate descrip- 
tion of many seemingly unrelated experimental findings, 
in terms of very few (theoretically calculable) parame- 
ters, namely: the k— mass mk [26] and the particle vi- 
bration coupling vertex h{abv), associated to the process 
in which a quasiparticle changes its state of motion from 
the unperturbed quasiparticle state a to 6, by absorbing 
or emitting a vibration v [1]. 

The equation describing the renormalization of a quasi- 
particle a, due to this variety of couplings is 
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where S ri and , (i ^ j) are the normal and abnormal 
self-energies. 

Eq. (l)is to be solved iteratively, and simultaneously 
for all the involved quasiparticle states. At each iteration 
step, the original quasiparticle states a with occupation 
numbers u a and v a , become fragmented over the different 
eigenstates a with probability u 2 a + v\, while the renor- 
malized occupation factors are obtained from the com- 
ponents of the eigenvectors, x a and y a , according to the 
relations u a = x a u a + y a v a , v a = -y a u a + x a v a . 
The quantities u a and v a are related to the spectroscopic 
factors measured in one-nucleon stripping and pick-up re- 
actions, respectively. One can also define [22,25] a renor- 
malized state-dependent pairing gap, through the rela- 
tion A a = 2E a u a v a /(u a + v 2 ), which in the limit of no 
fragmentation reduces to the usual BCS expression [27]. 

In the calculations reported below, a Skyrmc inter- 
action (Sly4 paramctrization, with mk ~ 0.7m[29]), 
was solely used to determine the properties of the bare 
single-particle states and the collective vibrations in the 
particle-hole channel. Both the bare nucleon-nucleon vu 
Argonne potential as well as the exchange of collective vi- 
brations were used in the particle-particle (pairing) chan- 
nel. 




FIG. 1: Renormalization processes arising from the particle 
vibration coupling phenomenon. A line indicates quasipar- 
ticles obtained from BCS theory, making use of the mean 
field single-particle states of Sly4 and the nucleon-nucleon v\4 
Argonne potential. The wavy line indicates the vibrational 
states. 



As seen from Fig. 3, Hartree-Fock theory is not able to 
account for the experimental quasiparticle energies of the 
low-lying states. Diagonalizing the Argonne vu nucleon- 
nucleon potential in the Hartree-Fock basis, within the 
framework of the generalized Bogoliubov-Valatin approx- 
imation including scattering states up to 800 MeV above 
the Fermi energy (to achieve convergence) in a spheri- 
cal box of radius equal to 15 fm, one obtains the state- 
dependent pairing gap shown in Fig. 4 (labelled Vu). 
The resulting pairing gap (average value for levels around 
the Fermi energy) accounts for about half of the empirical 
pairing gap value (« 1.4 MeV) obtained from the odd- 
even mass difference [30]. In keeping with this result, 
the quasiparticle spectrum (cf. Fig. 3), although be- 
ing slightly closer to the experimental findings than that 
predicted by Hartree-Fock theory, displays large discrep- 
ancies with observations. The situation is rather similar 
concerning the low- lying quadrupole vibration of 120 Sn 
calculated in the QRPA with standard effective nucleon- 
nucleon interactions like Gogny or Skyrme forces. While 
energy is predicted too high, which may not be too im- 
portant, the B(E2) value is about a factor 3 too small 
(cf. Table 1), a result which calls for a better theory. 

In fact, renormalizing the energy and the transition 
strength of the 2+ phonon, following NFT [2,3], that is, 
considering the couplings of the type depicted in Fig. 2 
(cf. also ref. [32] and refs. therein), one obtains an in- 
crease of the B(E2) transition probability which brings 
theory essentially in agreement with experiment (cf. Ta- 
ble 1) [33]. The most important processes which renor- 
malize the energy of the phonon are shown in Figs. 2(a) 
and (b). Other graphs which are also of fourth order in 
the particle-vibration coupling vertex, but contain inter- 
mediate states with more than four quasiparticle states, 
lead to very small contributions. This is because these 
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FIG. 2: Most relevant processes taken into account in the 
rcnormalization of the energy of the phonon (a-b) and of the 
associated transition strength (c-f). 
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FIG. 3: The spectra of the lowest quasiparticle states in Sn 
calculated using Hartree-Fock theory, BCS with the Argonne 
V14 potential, and after renormalization, are compared to the 
experimental levels in the odd neighbouring nuclei 119 Sn and 



terms not only involve larger denominators, but also, due 
to their higher degree of complexity, give rise to contri- 
butions with "random" signs which tend to cancel each 
other. This is a consequence of the fact that while can- 
cellation between the contribution associated wih graphs 
(a) and (b) of Fig. 2 is strong in the particle-hole channel 
the opposite is true in the particle-particle channel [35], 
and that the phonons are calculated in a Bogoliubov- 
Valatin-quasiparticle basis. In keeping with the above 
discussion, the most important processes renormalizing 
the B(E2) transition probability are those shown in Figs. 
2(c),(d),(e) and (f). 

We have also calculated the static quadrupole moment 
Q of the 2 + state, including the contributions from the 
processes shown in Fig. 6.27 of rcf. [1], considering also 
self-energy effects [36] . The resulting value of Q is rather 
small (8 e fm 2 ), in agreement with the experimental find- 
ings (10 ± 10 e fm 2 , or -5 ± 10 e fm 2 [37,38] ). 



FIG. 4: The state-dependent pairing gap for the levels close 
to the Fermi energy obtained using BCS theory with the «i4 
Argonne potential (circles) is compared with the result ob- 
tained including renormalization effects (squares) . 



Because in the above calculations we have included 
only a partial set (although the most important for the 
physics under discussion) of the NFT graphs needed to 
provide a completely consistent description of single- 
particle and collective vibration renormalizations, the 
mixing of spurious states with the physiscal states has 
to be contemplated. Although it is difficult to give a 
precise estimate of the error induced by such undesired 
couplings, 30% effects have been found in the calculation 
of the energy of the one-phonon state [39] . 

Making use of phonons which account for the exper- 
imental findings, the normal and abnormal self-energies 
were calculated, and Eq.(l) solved. The average value 
of the resulting state-dependent pairing gap of 120 Sn is 
now close to the value A. exp = 1.4 MeV derived from 
the odd-even mass difference (cf. Fig. 4) [40]. In Fig. 
3 we show the energy of the peaks carrying the largest 
quasiparticle strength, for the orbitals around the Fermi 
energy, which provide an overall account of the lowest 
quasiparticle states measured in the odd systems 119 Sn 
and 121 Sn. 

One can conclude that mean field theory and bare 
nucleon-nucleon potentials reproduce neither the exper- 
imental transition strengths nor the pairing gaps, let 
alone the density of quasiparticle states close to the 
ground state. Dressing the single-particle motion, the 
correlated particle-hole excitations of mean field and the 
nucleon-nucleon interaction with collective surface vibra- 
tions, brings theory in overall agreeement with experi- 
ment. In particular, about half of the pairing gap arises 
from the long range component of the pairing interaction 
associated with the exchange of collective vibrations. To 
further clarify the interdependence of single-particle and 
collective degrees of freedom, future studies should, for 
example, concentrate on the role this interdependence 
has on the nuclear masses. In particular, whether the 
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explicit, simplified, inclusion of ground state correlations 
and of the induced pairing interaction can reduce the 
present r.m.s. error of 0.674 MeV with which one of the 
best presently available Hartree-Fock mass formula [24] 
is able to reproduce the experimental findings. 
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hio 2+ (MeV) 


B(E2 T) (c" fm 4 ) 


RPA (Gogny) 


2.9 


660 


RPA (Sly4) 


1.5 


890 


RPA + renorm 


0.9 


2150 


Exp. 


1.2 


2030 



TABLE I: The energy and reduced E2 transition strength 
of the low-lying 2 + state, calculated according to different 
theoretical models, are compared to the experimental values 
[37]. 



